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Abstract 

Primordial fluctuations in inflationary cosmology acquire classical prop- 
erties through decoherence when their wavelengths become larger than 
the Hubble scale. Although decoherence is effective, it is not com- 
plete, so a significant part of primordial correlations remains up to 
the present moment. We address the issue of the pointer states which 
provide a classical basis for the fluctuations with respect to the in- 
fluence by an environment (other fields). Applying methods from the 
quantum theory of open systems (the Lindblad equation), we show 
that this basis is given by narrow Gaussians that approximate eigen- 
states of field amplitudes. We calculate both the von Neumann and 
linear entropy of the fluctuations. Their ratio to the maximal entropy 
per field mode defines a degree of partial decoherence in the entropy 
sense. We also determine the time of partial decoherence making the 
Wigner function positive everywhere which, for super-Hubble modes 
during inflation, is virtually independent of coupling to the environ- 
ment and is only slightly larger than the Hubble time. On the other 



hand, assuming a representative environment (a photon bath), the 
decoherence time for sub-Hubble modes is finite only if some real dis- 
sipation exists. 
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1 Introduction 



According to the inflationary scenario, all structure in the Universe originates 
from quantum vacuum fluctuations during a de Sitter (inflationary) stage in 
the very early Universe. These are inhomogeneous fluctuations (perturba- 
tions) of both the space-time metric and a scalar inflaton field. While the 
tensor part of the metric fluctuations generated during inflation produces the 
primordial gravitational wave background [lj, its scalar part together with 
perturbations of the inflaton leads to the origin of primordial density fluc- 
tuations producing present gravitationally bound objects and the large-scale 
structure of the Universe [2J. All these fluctuations can be formally repre- 
sented by a scalar field with a time- dependent 'mass', the time dependence 
coming from the coupling to the expanding universe described by a scale fac- 
tor a(t) (we assume here a spatially flat Friedmann- Robertson- Walker (FRW) 
cosmological model). 

Usually one assumes these fluctuations to be in their ground state (the 
adiabatic vacuum state) at the onset of inflation. This choice follows from 
the hypothesis of the maximal possible symmetry of the Universe in some 
period in the past during inflation. Also, it can be shown that this initial 
condition for fluctuations is an attractor for a wide open set of other initial 
conditions with a non-zero measure. The smallness of the fluctuations means 
that different modes (distinguished by a wave vector k with the wave number 
k = |k|) decouple and can be treated separately. Modes relevant for structure 
formation cross the Hubble radius H^ 1 twice. Here H = a/a is the Hubble 
parameter (with the dot denoting d/dt). During inflation, their wavelength 
A = 2na/k becomes bigger than H^ 1 ('the first Hubble radius crossing'). 
The quantum modes then acquire classical properties in the following sense. 
First, even without considering any interaction with other degrees of freedom, 
expectation values of any physical quantities constructed from the quantum 
fluctuations become practically indistinguishable from mean values of corre- 
sponding classical quantities as functions of classical stochastic fluctuations. 
This is achieved due to huge squeezing (in the sense of quantum optics) of 
those modes that gives a possibility to neglect non-commuting parts of all 
mode quantum operators. Since no consideration of environmental degrees 
of freedom and their interaction with modes involved are needed for this 
indistinguishability, this quantum-to-classical transition was called 'decoher- 
ence without decoherence' This result can be generalized to all higher 
non-linear orders of metric perturbations if only the growing (quasi-isotropic) 
mode of perturbations is kept jl]. In this way, it is possible to explain the 
classical stochastic behaviour of observed cosmological density fluctuations 
including their power spectrum and statistics. However, more subtle ques- 
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tions like the correct calculation of the entropy of primordial fluctuations 
cannot be solved in this approximation. 

Second, this classical behaviour is preserved, and even re-enforced, in the 
presence of an 'environment'. Since such an environment is always present in 
the form of other ('irrelevant') fields and fluctuations (or non-linear couplings 
of the same field), its influence has to be taken into account. In laboratory, it 
is usually the environment that leads to the emergence of classical behaviour, 
a process which is called decoherence (5j. In particular, strongly squeezed 
quantum states are known to be especially sensitive to decoherence. We 
have demonstrated some time ago that the classical basis distinguished by 
the environment (the 'pointer basis') is - in the limit of large squeezing - 
approximately given by the basis of field amplitudes [H] (see also [7J for a 
detailed conceptual discussion). The correct choice of the pointer basis is 
of particular relevance for the entropy of primordial fluctuations. We have 
calculated the entropy for each mode and found that it can assume at most 
half of its maximal value before the second Hubble radius crossing [8j . The 
maximal entropy would be achieved if the pointer basis was the particle- 
number basis [5]. 

The wavelength of the modes becomes again smaller than H~ l ('the sec- 
ond Hubble radius crossing') during the radiation and matter dominated 
phases. The scalar modes have left their imprint on the anisotropy spectrum 
of the cosmic background radiation, where they can be observed in the form 
of acoustic peaks (see e.g. PHE] for recent observations). Oscillations of the 
same origin (the 'Sakharov oscillations', often called baryon acoustic oscilla- 
tions) have been recently discovered in the three dimensional power spectrum 
of galaxy inhomogeneities [12] . These peaks would be absent if the entropy 
of each mode assumed its maximal value [5]. This demonstrates that ques- 
tions of the quantum-to-classical transition for primordial perturbations are 
not only of academic nature, but have observational relevance. Primordial 
gravitational waves (tensor modes) generated during inflation also produce 
primordial peaks in the CMB temperature anisotropy and polarization mul- 
tipoles (see [13] for discussion of perspectives of their observation). 

In spite of these investigations, the discussion about the precise mech- 
anism of this quantum-to-classical transition and the amount of entropy of 
primordial fluctuations gained in its course goes on, cf. the recent papers 
[El EEi EES CEO HH1 EH]. It was argued, for example, that the pointer basis 
is not given by the field-amplitude basis, but by the basis of two-mode co- 
herent states [H] , by the squeezed-state basis [IB], or by the second-order 
adiabatic basis [20J . It was further argued in [H] that the associated entropy 
of fluctuations (being half the maximal entropy) provides a lower (instead of 
upper) bound on the entropy. Also, there is a disagreement between [18] and 
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[12] regarding the question whether decoherence may occur already during 
inflation (after the first Hubble radius crossing, of course) or only after its 
end, during reheating. We find it therefore appropriate to address this issue 
again by reviewing and extending our line of arguments. 

Our article is organized as follows. In Section 2 we first state the problem 
of the quantum-to-classical transition and briefly review the central general 
concepts. We then apply these concepts to the primordial fluctuations. We 
give exact and approximate expressions for the von Neumann entropy and 
also calculate the linear entropy. We then discuss why the pointer states for 
the fluctuations is given by approximate field-amplitude eigenstates (narrow 
Gaussians) and not by another basis such as the coherent-state basis. In 
Section 3 we apply methods from open-system quantum theory (the Lind- 
blad equation) to calculate in a general setting decoherence times for modes 
outside and inside the Hubble radius. We conclude from this again that 
the pointer states are narrow Gaussians in the field amplitude. The entropy 
per mode is thus bounded from above by half the maximal entropy in the 
long-wave regime before the second Hubble radius crossing. 

2 Quantum-to-classical transition and entropy 

2.1 Quantum description for the primordial fluctua- 
tions 

The primordial fluctuations which lead to the observed structure in the Uni- 
verse can originate from the quantum fluctuations of the metric and a scalar 
field in an early phase where the Universe is expanding in a (quasi-) expo- 
nential way ('inflation'). These fluctuations are assumed to be small so that 
their mutual interaction can be neglected in the description of the main sce- 
nario. This linear treatment of the fluctuations is certainly justified by the 
smallness of the observed CMB anisotropy and of the ratio of the modulus 
of the Newtonian potential to the square of the light velocity for all observed 
gravitationally bound structures in the Universe (apart from the very close 
vicinity of black holes). It is thus appropriate to represent the modes in 
Fourier spaces where the various modes (wave numbers) decouple form each 
other. The dynamics of both the tensor and the scalar modes can be repre- 
sented by some scalar field 0k- Because of the mode decoupling we shall in 
the following skip the index k. It turns out to be convenient to work with 
the rescaled field amplitude y = a<p. In the case of tensor fluctuations (grav- 
itational waves), this is already the appropriate variable to deal with; for 
scalar perturbations one has to use a gauge-invariant combination of metric 
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and inflaton perturbations. These details are irrelevant for the study of the 
quantum-to-classical transition. 

Since we are working in Fourier space, y is complex. This reflects the 
fact that we shall have a two-modes state, the two modes being given by 
k and — k. For the following discussion it is, however, sufficient to assume 
that y represents one real mode, but one should keep in mind that the actual 
number of degrees of freedom is twice as much. This also holds for the entropy 
discussed below; the entropy for the two-mode state is twice as much as for 
the one mode considered here. 

The dynamics is then described by the Hamilton operator 

h = \(p 2 + *Y + ^«p) , (i) 

where a' = da/drj, with rj denoting conformal time, and p is the momentum 
canonically conjugate to y. (We have skipped the index k.) This Hamiltonian 
follows after expanding the metric and the scalar field up to second order 
in the inhomogeneities, with the cosmological background given by a flat 
Friedmann-Lemaitre universe. 

Each mode satisfies the Schrodinger equation 

«^=A«. (2) 

Assuming that the modes are initially (at the onset of inflation) in their 
ground state, their wave function is then at any time of the Gaussian form, 

VKv. v) = (^v^) V4 exp (-my 2 ) , n = n R + in, . (3) 

Non-Gaussian initial states can also be studied [21]; it was shown there that 
there is a wide class of initial non- vacuum states that cannot be distinguished 
from a vacuum initial state by just looking at the statistics of observable 
quantities. 

One can prove from the form of the Hamiltonian (TjQ) that the states (EJ) 
represent squeezed states; the squeezing is generated by the last term in (DO). 
(Actually, we have a two-modes squeezed state, but for simplicity we present 
one mode only, cf. the remark above.) Introducing the squeezing parameter 
r and the squeezing angle ip, one can write (setting h = 1 from now on) 

k 

tt R = — — ■ . — , fli = -fi R sin2y2sinh2r , (4) 

cosh 2r + cos 2(p smh 2r 
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with = k and fii = for the initial vacuum state. Both r and (p are 
functions of t (resp. 77), the exact dependence arising from the expansion law 
a(t). For pure exponential inflation described by 

i 

a(t) = a exp(Hjt) - 



( V - 27 7e )# I ' 

where r] e denotes the conformal time at the end of inflation, simple analytic 
results are available corresponding to the adiabatic vacuum in the de Sitter 
space-time [31 [22] : 

sinhr = — — , cos2<£> = tanhr . (5) 

This would lead to r — > oo, tp — > if inflation lasted arbitrarily long. Infla- 
tion, however, ends, but it still leads to r « 120 for the largest cosmological 
states [23], so one can certainly neglect after some time terms of the order of 
exp(— r) ('decaying mode'). From (Hj) one then finds the asymptotic values 

tt R -> ke~ 2r , Qi -> -ke~ r , (6) 

exhibiting f2j 3> (WKB limit). This behaviour is assumed to occur gener- 
ically for an inflationary model. One recognizes that the width of the Gaus- 
sian ([3]) becomes very broad in y and highly squeezed in p. This broadness 
reflects the fact that the kinetic term of the Hamiltonian becomes negligible 
in this limit; in the Heisenberg picture this is exhibited by the fact that the ki- 
netic term of the Hamiltonian, H^i n , then approximately commutes with the 
field-amplitude operator, y. In the limit of high squeezing, quantum expecta- 
tion values become indistinguishable from mean values of classical stochastic 
quantities, the difference containing typically terms of the order exp(— r) or 
powers thereof. In a quasi de Sitter stage, the annihilation operator evolves 
according to 

a(k)-^(a(k) -at(-k) ) (7) 

for k -C aH. So far, the primordial fluctuations behave fully quantum. They 
are described by states which are very broad wave packets. Such broad 
packets would exhibit in a laboratory situation very non-classical behaviour, 
as could be checked in interference experiments. The standard scenario of 
structure formation starts, however, from classical stochastic fluctuations. A 
major issue is thus to describe the quantum-to-classical transition for the 
primordial modes. We shall now first briefly review the general concepts 
and subsequently shall apply them to understand the emergence of classical 
behaviour for these modes. 
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2.2 General concepts of the quantum-to-classical tran- 
sition 

The superposition principle is at the heart of quantum theory. As a conse- 
quence, the set of classical states is of measure zero, since one can always 
superpose different classical states (e.g. different narrow wave packets) to 
get 'weird' non-classical states. 

Such states may occur in many ordinary measurement situations. If we 
assume that a measured system is initially in the state |n) and the mea- 
surement device in some initial state |$o) ; the evolution according to the 
Schrodinger equation in the simplest case reads 

|n)|$ ) exp (-iH int t) |n)|$ > = |n)|$„(*)) , (8) 

where H mt is a special interaction Hamiltonian (assumed here to dominate 
over the free Hamiltonians) which correlates the system state with the device 
without changing the former. The resulting apparatus states |$ n (t)) are often 
called 'pointer states'. In order to yield distinguishable outcomes, these 
pointer states must be approximately orthogonal. A process analogous to 
(JHJ) can also be formulated in classical physics. The essential new quantum 
features now come into play when one considers a superposition of different 
eigenstates (of the measured 'observable') as the initial state. The linearity 
of time evolution immediately leads to 

(j2c n \n) J |$ ) J2cn\n)\$ n (t)) . (9) 

\ n / n 

But this state is a superposition of macroscopic measurement results. It 
is by now well established, both theoretically and experimentally, that the 
ubiquitous and unavoidable interaction with the environment has to be taken 
into account |5J. The measurement device is itself 'measured' (passively 
recognized) by the environment according to 

(^c»|$ n > J \E ) ^c»|$ n )|£ n >. (10) 

\ n / n 

This is again a macroscopic superposition, now including the myriads of 
degrees of freedom pertaining to the environment (gas molecules, photons, 
etc.). However, most of these environmental degrees of freedom are inac- 
cessible. Therefore, they have to be integrated out from the full state ( iTUl) . 
This leads to the reduced density matrix for system plus apparatus, which 
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contains all the information that is available there. It reads 

Psa ~ Yl M 2 | n )H ® l $ n)( $ ™| if (E n \E m ) S3 5 nm , (11) 

n 

since under realistic conditions, different environmental states are orthogonal 
to each other (they can discriminate between different states of the appara- 
tus). Equation fill I) is identical to the density matrix of an ensemble of 
measurement results |n)|$„). System and apparatus thus seem to be in one 
of the states \n) and |$ n ), given by the probability |c„| 2 . This is decoher- 
ence. The interaction with the environment singles out a preferred basis for 
the apparatus, which is just the above-mentioned pointer basis. 

Assuming the usual Born rule for quantum probabilities, the reduced 
density matrix contains all the information that can be obtained from system 
and apparatus alone, without taking into account the quantum entanglement 
with the environment explicitly. 

In our case it is the primordial fluctuations for which a classical behaviour 
has to be obtained. Formally, they correspond to the apparatus states above, 
but since we have no additional system, we just call them the system vari- 
ables; they are in interaction with the environment. The latter is represented 
by other fields or by higher-order modes of the metric and inflaton itself. 

An important property of a pointer basis is its preservation in time - they 
should be robust while interacting with the environment, that is, they should 
not evolve into superpositions of themselves. This basis thus represents the 
properties that become classical. For this, two conditions have to be fulfilled 
[5]: first, the projection operators on the pointer states should approximately 
commute with the interaction Hamiltonian between system and environment; 
second, these projection operators should approximately commute at differ- 
ent times (in quantum optics, this is referred to as a quantum non-demolition 
condition). We shall see in the following subsections that these conditions 
are fulfilled in the high- squeezing limit for the field-amplitude basis of the 
primordial fluctuations. 

In the literature, various methods are discussed with which one can deter- 
mine the pointer basis [3] . It has been shown that they are equivalent to each 
other at least in simple situations [24J . One prominent method makes use of 
the entropy of the system. As long as it is in a pure state, this entropy is 
zero. If it is in a mixed state due to decoherence, this entropy is positive. Its 
value characterizes the degree of entanglement between system and environ- 
ment - the stronger the entanglement, the higher the entropy. One can thus 
find the most robust system states (that is, the pointer states) by minimizing 
the resulting entanglement entropy. This method is called the 'predictability 
sieve' [25] . We shall apply this method below to the primordial fluctuations. 



7 



A final tool for the study of the quantum-to-classical transition should 
be introduced - the Wigner function. For a general density matrix in the 
position representation, p(x,x'), it is defined by the expression 



W(x,p) = — / dy e 2tpv p(x - y,x + y) 

7T / 



1 f°° 



(12) 



J — oo 



It is thus defined on phase space where it gives the correlations between 
position and momentum. It is not, in general, positive definite, reflecting the 
fact that quantum theory is inequivalent to a classical statistical theory on 
phase space. 

2.3 Wigner function and reduced density matrix for 
the primordial fluctuations 

We shall first address the system of primordial fluctuations as if it were 
isolated [3]. In order to recognize correlations between position (here: field 
amplitude) and momentum (here: the variable canonically conjugate to the 
field amplitude), the Wigner function can be calculated. In general it can 
assume negative values, but for the special case of a Gaussian wave function 
it is everywhere positive (of a Gaussian form) . This is the case here where the 
fluctuations are assumed to be in their ground state. A convenient quantity 
is then the 'Wigner ellipse' defined by the contour of this Gaussian Wigner 
function, because it provides one directly with a measure for the mean square 
deviations of y and p. From the analysis of the Wigner function one can see 
that the system behaves as if it had a stochastic amplitude (with a Gaussian 
distribution if the state is (j3J)) and a fixed phase, see the calculations in 
[3J El [7]. This phase is given by the squeezing angle. The situation here is 
thus already a peculiar one since the quantum nature of the squeezed state 
cannot be seen, in the high-squeezing limit, from expectation values and 
variances [3]. 

The scenario for the isolated system of primordial fluctuations proceeds 
as follows. After the first horizon crossing (during inflation) one has r < 100, 
(p — > 0, and the Wigner ellipse will have - after an appropriate re-scaling of 
the axes - a major half axis a = exp(r) and a minor half axis f3 = exp(— r). 
After the second horizon crossing (in the radiation or matter dominated 
phase) the Wigner ellipse rotates, with small oscillations around a large mean 
value of r. The rotation is very slow, the frequency being about the inverse 
age of the universe, which is why the squeezed nature is preserved for a long 
time [21110 For modes that re-enter in the radiation era, this rotation leads to 




In contrast to this, the rotation is fast in the case of black holes that has important 
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the acoustic peaks in the cosmic background radiation [31 [22]- These peaks 
have been observed with high precision [TUJ [TT] . 

This discussion treats fluctuations as being isolated, which is, however, 
unrealistic [HIEEJEH]. Other fields (the 'environment') interact with them, 
even if the coupling is weak. Such a coupling may also arise from a small 
self-interaction of the modes [T3]. Even if there were no other fields present, 
this small self-interaction would remain. As long as the modes are outside the 
Hubble radius, this coupling cannot lead to a direct causal interaction, but 
can only produce entanglement ('EPR-type situation'). Since usual interac- 
tions couple to fields (instead of canonical momentum), the coupling is in y 
(not p). This already suggests that the y-basis is equal to the pointer basis 
at least approximately. Spatial gradients would not change this conclusion, 
since they do not depend on p. 

The situation is analogous to the localization of particles in quantum me- 
chanics (5j [30]. This interaction with the environment without direct causal 
contact can be described by multiplication of the density matrix p (y,y') 
corresponding to the system alone with a Gaussian factor according to 

po(y,y') — >pdv,y') =Po(y,y') ex P (-^(y-y') 2 j > ( 13 ) 

cf. 0, ESJ. The parameter £ encodes phenomenologically the details of the 
interaction strength with the environment. This interaction is by no means 
restricted to be linear and can be very complicated. In this way, non-diagonal 
elements are suppressed with respect to the y-basis. In a realistic situation 
for decoherence, one would expect that £ dominates over the corresponding 
part in ©, 

£ » n R « ke~ 2r . (14) 

The typical time scale for decoherence during inflation is ~ H^. (There 
is a close analogy to chaotic systems, with Hi corresponding to the Lyapunov 
parameter [8]. The reason for this is the classical instability of our system.) 

From analogous investigations in quantum mechanics [31] . it can realis- 
tically be expected that the parameter £ settles to a constant value after a 
certain transition time. The axes of the Wigner ellipse then read 




where ( JT41) has been used. While the major axis has remained unchanged, 
the minor axis has become bigger and settles to a constant value. The area 

consequences for Hawking radiation [27] . 
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of the Wigner ellipse thus increases, leading to a non-vanishing entropy (see 
below). In order to avoid conflict with observation, one has to impose the 
'correlation condition' (3 <C a, leading to 

f < e 2r . (16) 
k 

Together with the slow rotation of the ellipse, this guarantees the formation 
of acoustic oscillations. A random distribution of p and q would totally smear 
out these structures. 



2.4 Entropy 

As long as the fluctuations are treated as an isolated system, they can be 
described by a wave function and thus have vanishing entropy. This is no 
longer the case for the interacting system described by ( fl3i) ; the entanglement 
with the environment leads to a 'loss of information' for the system and 
therefore to a positive entropy. The quantum correlations are present only in 
the total system and are therefore unseen 'locally'. The entropy is calculated 
from 

S = -trfehxpe) , (17) 

where we have set &b — 1- The entropy has, of course, also to be calculated 
for a two-mode squeezed state, although for simplicity we give here again the 
calculation for one real mode only. For the full entropy per the two-mode 
state, one thus has to double the results below. 

It is convenient to introduce the dimensionless parameter \ = £,/Qr con- 
trolling the strength of decoherence. (In the case of pure exponential inflation 
one has x — £(1 + 4 sinh 2 r)/k.) Inserting (Tl3|) into ffl7|) . one gets the explicit 
expression 



s = -in =L ifv^-iWV 1 ^- 1 



+ X + 1 2V V ~ J v/T + ^+1 

ln oVX- V 1 + Xln — ■ 

2 y/X 



18) 



One recognizes that the entropy vanishes for £ — > 0, as it must. In the limit 
X 3> 1 (large decoherence) one gets 

S=l-ln2 + ^ + 0(;T 1/2 ) • (19) 
Both ( {TBI) and ( 1191) are displayed in Figure 1. 
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Figure 1: Entropy in dependence on the decoherence parameter x- Shown 
are the exact expression and the limiting case. 

It is seen that the asymptotic value is readily attained. For \ ~ 1-5 
one obtains S ~ In 2, corresponding to the loss of one bit of information. 
Usually decoherence sets in roughly at this stage [5j. A stronger condition 
than (TIE]) is obtained by the reasonable demand that some squeezing should 
remain compared to the vacuum state which has Or = k. This leads to the 
condition £ < k and yields for exponential inflation in the high-squeezing 
limit the bound S < r. Whether this really holds is, of course, a question 
about realistic interactions in the universe. 

The maximal entropy is instead given by S max = 2r (obtained by smearing 
out the Wigner ellipse into a big circle and corresponding to the pointer basis 
being the particle- number basis [§]). One can thus define a notion of partial 
decoherence in the entropy sense by the ratio of entropy to maximal entropy, 
S/Smax., P er eac h mode. 

It is also instructive to consider the linear entropy. It is defined by 

Siin = tr (p - p 2 ) , (20) 
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obeying < S\ in < 1. Inserting ([TBI , one obtains 

s lin = i - (i + xr 1/2 • (21) 

For x ^ 1 one gets SW ~ x/2, while for x ^> 1 one gets Sy m wl - X~ 1//2 < 
1 — exp(— r), where we have used the bound £ < k in the last step. 

As was already mentioned, the fact that the interaction with the envi- 
ronment is in y-space leads to the suggestion that - for large squeezing - the 
pointer basis equals approximately the field-amplitude basis [6]. Exact diago- 
nalization would correspond to an entropy increase given by S = r = S max /2 
[32"] . in accordance with the expectation that S < r holds for modes outside 
the horizon. Contrary to this it was claimed in [14J that the pointer basis 
is given by (two-mode) coherent states. Diagonalization in such a basis was 
in the cosmological context first addressed in [33], with the result that S 
approaches for large squeezing the value r = S max /2, that is, the same value 
as for the field-amplitude basis. Figure 2 displays the entropy S(r) for the 
field-amplitude basis, the coherent-state basis ('^-basis'), and the particle- 
number basis ('n-basis'). It is easily seen that S approaches r in the first two 
cases, whereas it approaches 2r in the latter case. 

2.5 Pointer basis for the primordial fluctuations 

We shall now re-enforce from various points of view our earlier result P, [7] 
that the pointer basis is for large squeezing (approximately) given by the 
field-amplitude basis. 

The first argument is about the normalizability of the pointer basis. The 
authors in [H] are concerned about the fact that field-amplitude eigenstates 
are not normalizable and that they, moreover, exhibit infinite momentum 
spread. However, even if one dealt with exact amplitude ('position') eigen- 
states, the quantum mechanical formalism would be able to cope with such a 
situation (GNS construction instead of Hilbert space). But the point is that 
the environmental parameter £ entering (1131) is never strictly infinite, so that 
a finite, albeit small, off-diagonal part for the density matrix remains. One is 
thus dealing with very narrow wave packets that are perfectly normalizable. 
Even if we had exact position eigenstates, an 'infinite' momentum spread 
would only take place if the kinetic term played a crucial role in the dynam- 
ics. As has already been remarked above, however, for the highly squeezed 
modes the kinetic term is negligible compared to the potential term. 

This irrelevance of the kinetic term is also the reason why the reference 
to [25] in [2] is misleading. In [25], the question of the pointer basis has 
been studied for an ordinary harmonic oscillator - in contrast to cosmology 
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Figure 2: Entropy S in dependence of the squeezing parameter r for three 
cases of the pointer basis. 
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where one deals more with an upside down oscillator [33] - coupled to a 
high-temperature environment, using the method of the 'predictability sieve' 
(minimal local entropy production). It was found that this basis is given 
by the coherent-state basisU A crucial ingredient in this proof, however, 
was the importance of both the kinetic and potential term for the oscillator 
(an averaging procedure over many oscillation cycles had to be performed). 
For small times, when the kinetic term does not yet become relevant, it was 
found in [25] that, in fact, the position basis is the (approximate) pointer 
basis. This is the limit of relevance for the cosmological fluctuations. 

In situations where the system Hamiltonian is negligible, the pointer basis 
has the property that it consists of states which are eigenstates of an operator 
that commutes with the interaction Hamiltonian. If the system Hamiltonian 
is not negligible (in our case, this happens for small r), one must invoke a 
principle such as the 'predictability sieve' or the 'rate of de-separation' in 
order to determine a pointer basis. In the case of a free particle coupled to a 
localizing environment it was shown that the predictability sieve (as well as 
other criteria) predict robust pointer states that are narrow Gaussian wave 
packets [24] . In the next section we shall apply methods of open system 
dynamics to the Hamiltonian (CQ) in order to discuss the decoherence time 
for the cosmological fluctuations with wavelengths both bigger and smaller 
than the Hubble scale. This will enforce our earlier result about the nature 
of the pointer basis. 

Another argument follows from the requirement that the pointer basis 
should be stable with respect to time evolution. In the limit of no interac- 
tion with environment, this means that it should be stable with respect to 
free evolution. However, the coherent state basis does not satisfy this re- 
quirement in the period between the two Hubble radius crossings: expansion 
of the Universe causes coherent states to become strongly squeezed ones in 
a characteristic time ~ H . 

This can be seen as follows. It has been shown [221 136] that the time 
evolution generated by any quadratic Hamiltonian like (OQ) factorizes into an 
irrelevant phase rotation and the action of a squeezing operator with param- 
eters r and <p. It is further well known (e.g. [37J Sect. 2.7, and [36J ) that a 
squeezing operator acting on an arbitrary coherent state \z) (not just the vac- 
uum ground state) produces a 'squeezed coherent state': It has characteristic 
uncertainties proportional to e~ r and e r along rotated axes in the y-p space, 
coinciding with the squeezed and the stretched semi-axes of the Wigner el- 
lipse, respectively. If we thus consider the time evolution of a coherent state 

2 This result had already been found by another method - implementing the rate with 
which initially unentangled states become entangled - in |35j . 
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starting at a time t\ after the mode has left the Hubble scale, we have to be 
concerned solely with the development of the squeezing parameters r and (p 
starting with r = = ip at that time. As in the particular de Sitter example 
mentioned before, the squeezing angle <p then quickly 'freezes out' and the 
squeezing factor grows like r pa ln(a/aj), a result first obtained for gravi- 
tational waves [23J, and later generalized to scalar perturbations [22] • For 
quasi-exponential inflation, we consequently have r pa Hi(t — ti), and it takes 
a very short time ~ iff 1 for the state to become significantly squeezed. But 
even if we start in the ensuing radiation-dominated era, where a oc {t/t e ) l l 2 
and H = l/(2t), we obtain r pa m(t/ti)/2 = \n(2Hit)/2, so that squeezing 
still becomes effective at a time t ~ (e 2 /2)H^ 1 = e 2 ti (with H\ denoting the 
Hubble parameter at the initial time tC). Finally, these considerations also 
apply to two-mode coherent states (cf. the above references). On the other 
hand, the field-amplitude basis is stable with respect to time evolution in the 
super-Hubble regime, since the field amplitude remains constant with great 
accuracy. Here, the squeezing manifests itself in the continuous decrease of 
the decaying mode (see e.g. [35]). 

An interesting and subtle situation arises with the second order adiabatic 
basis proposed in [20]. In the case of the exactly massless and minimally cou- 
pled scalar field (CQ) and the exact de Sitter background, the vacuum mode 
functions (f) k of the second order adiabatic vacuum (as defined in this paper) 
coincide with the exact solution for <ft k = y k /a. This has led the authors of 
[2"U] to the statement that there is no creation of massless minimally cou- 
pled particles in the de Sitter space-time. However, this result is a purely 
academic one since in case of the stable, eternal de Sitter expansion there 
is no possibility for an observer to measure these particles because they are 
beyond his/her future event horizonH On the other hand, for a viable infla- 
tionary model in which the Hubble parameter H decreases during inflation, 
the mode functions <p k do not coincide with the exact solution for <p k which is 
constant far outside the Hubble radius (k <C aH) and equal to H{t k ) / V2k 3 
there (t k is the moment of the first Hubble radius crossing when k = aH). 

Really, the second-order adiabatic vacuum basis functions were defined 
in [20] in the form 

exp l-i / dt'W k {t') \ , Wk = -— . — (22) 



y/2a?W k V. J J a 2k 2 + a 2 (H + 2H 2 ) 

for a spatially flat FRW model with an arbitarty a(t). For a quasi-de Sitter 
(slow-roll) regime \H\ -C H 2 , a(i]) pa —1/H(i])r] (77 < 0), this expression 



3 So, it is reminiscent of the problem, much discussed in the past, of radiation reaction 
for a charged particle eternally moving in a constant electric field. 
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reduces to 

Outside the Hubble radius, \<j) k \ = H(r])/VW oc H(t). Since H(t) < H{t k ) 
for t > ifc, this means that creation of massless minimally coupled particles 
during slow-roll inflation does occur even with respect to this basis. By the 
end of inflation and subsequent heating of matter, we return to standard 
formulas for the particle number and the energy density power spectrum 
presented in [T] (for the case of gravitons). Thus, in contrast to the field- 
amplitude basis, the second-order adiabatic basis (1221) is also not stable with 
respect to time evolution in the super-Hubble regime even during an infla- 
tionary stage with slow-roll. 

Of course, there exists an ambiguity in the definition of adiabatic vacuum 
mode functions of a given order (in other words, in the definition of the notion 
of number of particles in an external varying gravitational field) - alternative 
definitions may differ by terms of a higher order in the adiabatic parameter 
aH I k and its time derivatives. However, different variants of adiabatic mode 
expansion produce essentially the same form of the adiabatic expansion of 
the average value of the quantum field energy-momentum tensor as was 
shown already in [39l HQ] . Note in this connection that the most standard 
textbook calculation of a next small correction to the WKB solution of the 
wave equation 

k 2 

<P k + 3H(t)(f> k + —--<j) k = (24) 
a z (t) 

(that corresponds to the Hamiltonian (p"}) in the WKB (sub- horizon) regime 
k aH) reads 

fe = * e ~*7(l + /) , f = 7 Lf dt a(H + 2H 2 ) 1 (25) 
ayzk J 

If we formally use (125}) for all k during slow-roll inflation ignoring the condi- 
tion |/| -C 1, it just reduces to the same expression ([23]) . exact for H = const 
but significantly different from the exact solution for a slowly variable H. See 
also (HI H21 S3] for recent further improvements of the WKB-expansion in 
this case, in particular, making it homogeneous in both the sub- and super- 
Hubble regimes. 

We would finally like to address the general criticism of |44j where it is 
claimed that in the discussion of emerging classical behaviour for the fluc- 
tuations a transition from a symmetric quantum state to a non-symmetric 
classical state is implicitly assumed without justification. The authors of 
[4*4"] therefore claim that in addition new physics (an explicit wave function 
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collapse) is needed. This is, however, a general issue in quantum theory 
and plays a role, for example, in spontaneous symmetry breaking, see [45J, 
Sec. 6.1: the initial symmetric state develops into a symmetric superposition 
of all 'false vacua'. The standard false vacuum is obtained by selecting one 
component out of this superposition that can be justified using any interpre- 
tation of quantum mechanics, for example the Everett or Copenhagen one, 
without changing of quantitative predictions of quantum mechanics referring 
to this component. The same happens here: the initial symmetric vacuum 
state evolves into a symmetric superposition of inhomogeneous states out of 
which one component is 'selected' [35]. Thus, cosmological perturbations are 
not specific in this sense and neither solve nor complicate the fundamental 
problem of the foundations of quantum mechanics. 

3 Decoherence time and pointer states 

3.1 Description of decoherence by master equations 

Decoherence is the irreversible emergence of classical properties for a quan- 
tum system through its unavoidable interaction with the 'environment', that 
is, with irrelevant degrees of freedom [5] . Ideally, one would solve the Schro- 
dinger equation for system plus environment and then trace out the envi- 
ronmental degrees of freedom to obtain the reduced density matrix for the 
system. This density matrix obeys a master equation which provides a non- 
unitary and irreversible time evolution. Instead of performing this procedure 
explicitly, one can directly make a general ansatz for the master equation 
that can cope with all interesting situations. A convenient form for this 
equation is the 'Lindblad form', cf. Sect. 3.3.2.2 in [5] or Eq. (7) in [36]. The 
corresponding master equation is Markovian (local in time) and preserves 
the properties of a density matrix (such as conservation of its trace). Such 
an equation results from a wide range of realistic interactions; the details 
of these interactions are encoded in the Lindblad operators. The simplest 
form is a pure localization term in addition to the system Hamiltonian; such 
a term leads to the Gaussian suppression of interferences described by (TT3l . 
The most general interaction may be non-Markovian, but since we are inter- 
ested in the minimal mechanism to guarantee decoherence, the restriction to 
the Markovian case is sufficient. 

In [31], a master equation for the reduced density operator p was studied 
for a free particle plus such a localization term. It reads 

^-ass^-ffcfcwn- (26) 
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Here, D contains the strength of the interaction with the environment and 
describes the strength of localization. The term in f[2"oT) containing D can 
arise from a high-temperature environment, as it was discussed for example 
in [25], but it does not have to. It actually follows from a much wider class 
of situations [30]. Even if the environment is thermal, the temperature does 
not need to be high; in a well-known example calculated in [30] a small 
dust particle in intergalactic space is localized (decohered) by the ubiquitous 
microwave background radiation, a process described by an equation of the 
form (1261) . This master equation is also obeyed by the density matrix for the 
primordial fluctuations, which is discussed in [8]. 

We shall here consider the situation where a general initial state (not 
necessarily Gaussian) is present. In such a situation the Wigner function 
is usually not positive definite (cf. [IB] for the realization of such states 
in the laboratory). The emergence of positivity is usually connected with 
decoherence, so it can be used alternatively to the density matrix as a measure 
for classicality. It was found in [I7J that the Wigner function for the case 
( T2T)]) becomes positive after a certain time t^, independent of the initial state. 
This, therefore, signals the emergence of classical behaviour. Moreover, it 
was shown there that the reduced density operator p can for t > t& be 
decomposed in the form, 



where |T) denotes a set of Gaussian states. They play the role of the pointer 
states, so ( 1271) is not the standard orthogonalization of the density matrix; 
note that the decomposition ( 1271) is with respect to an overcomplete basis. 
Localization leads here to the emergence of narrow wave packets in position 
space (but not to delta functions). The width of these narrow Gaussians was 
determined by the 'predictability sieve' [21]. Apart from numerical factors 
of order one, the width is given by the expression 5 = (h/mD) 1 ^. If the 
kinetic term in (1261) is neglected compared to the localization term, we have 
m — > oo and thus 5 —>■ - the width of the Gaussian pointer states becomes 
very small, but they remain of course normalizable. For a small kinetic term 
the pointer states are thus approximate position eigenstates. This is the 
situation that has its analogue in the long-wavelength modes in cosmology. 
The limit m — > oo is analogous to the limit r — > oo. This can be seen from 
the correlation p « (tan</?)y « e~ r y, which holds for exponential inflation 
and expresses the smallness of the kinetic (p 2 -containing) term. 

We shall address below the question of positivity for the Wigner function 
and the ensuing decoherence time in the cosmological case. While for the 
special initial state ([3]) the Wigner function is always positive, the question 




(27) 
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of positivity is non-trivial for all other initial states, such as the ones in [2"T| 
(coherent states, states with an arbitrary number N of particles). 

The results of [U] were generalized in [16] to a system Hamiltonian of 
arbitrary quadratic form and to a general Lindblad equation. It was shown 
in particular that the Wigner function has always (apart from the case of 
Gaussian states) negative parts for t < t d . The emergence of positivity at 
time td (and not earlier) is thus independent of the initial state. It was 
generalized to non-Markovian situations in [19]. It can thus be assumed that 
in a generic situation of a particle coupled to an environment, the state of 
the system can after a finite time not be distinguished from an exact mixture 
of particular Gaussian states which have a narrow width in case of strong 
interaction. 

3.2 Lindblad equation for the primordial fluctuations 

We shall restrict ourselves in the following to the case of one Lindblad oper- 
ator L. The generalization to several Lindblad operators is straightforward 
[S] . Such a generalization is needed if one wants to accommodate non-linear 
effects by the primordial modes themselves, that is, to express the influence 
of 'environmental' fc-modes on 'system' k- modes. Such a master equation 
was used in [TBI 02] , whereby the 'system modes' were chosen to be the long- 
wavelength (observable) modes and the 'environmental modes' were chosen 
to be the short-wavelength (unobservable) modes with wavelengths smaller 
or of the order of the Hubble radius. 
The master equation then reads [5] 

^ = -i[H, p\ + LpD - \tihp - l -pDL . (28) 

We further assume as in |46j that the Lindblad operator is linear in our 
variables p and y. The general system Hamiltonian discussed in [IB] is of the 
form 

H = H u p 2 + 2H 12 py + H 22 y 2 . (29) 

To conform with the conventions used there, we will from now on work with 
a slightly different set of canonical variables, defined by the substitutions 
p — > \fkp and y — > y/ \f~k~ in (JTJ), so that the Hamiltonian reads 

H = ^{ P 2 + y 2 ) + %p. (30) 

In our case we thus have H\\ = H 22 = k/2, and H\ 2 = a' /2a = T~i/2, where 
Ti = aH denotes the conformal Hubble parameter. In the analysis of [46J, the 



19 



determinant of the quadratic form defined by (|2"5]1 plays a crucial role. (We 
note that this analysis relies on a Hamiltonian that does not explicitly depend 
on time. Strictly speaking, therefore, our treatment is only valid for times 
that are reasonably short compared to the characteristic timescale given by 
TC.) More precisely, the sign of this determinant decides about the qualitative 
features of the decoherence time. In our case, this determinant is positive 
for k > aH = TC ('elliptic case') and negative for k < aH ('hyperbolic 
case'), that is, the change of sign just happens when the modes cross the 
Hubble radius. For modes outside the Hubble radius, the hyperbolic case 
is of relevance. This turns out to be of crucial importance. The 'inverted 
oscillator' is the paradigmatic example for the hyperbolic case |34j . 

Further following the notation used in [IB] , we write the Lindblad operator 
in the form 

L = (l' + il").Q, (31) 
where the real components of the 'vectors' 

r= © and '" = (y) (32) 

contain all information on the interaction of the system with its environment. 
We also introduce the 'dissipation coefficient' a = //A" — //'A' as well as the 
quantity a = 2y/- det H = y/TC 2 - k 2 . Note that cr is real for k < TC (outside 
the Hubble radius) and imaginary for k > TC (inside the Hubble radius) . The 
parameter a is non-vanishing if there is an energy exchange between system 
and environment. The real part of a is called 'Lyapunov exponent' in [46J, 
which is in accordance with our discussion of the entropy in [8]. It is non- 
vanishing only in the hyperbolic case. The quantum-mechanical example 
(126p is recovered from ( 128]) in the special case of A' = A" = 0, and one can 
then identify D = fi' 2 + fi" 2 . This leads in particular to vanishing friction, 
a = 0. 

A quantitative result of following |47j, is that the positivity time, i] p , 
for the Wigner function is the solution of 

detM(-»j) = i, (33) 

with a certain (time-dependent) matrix M that appears in the solution for 
the Wigner function and that depends on the Lindblad operators. For the 
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case of primordial fluctuations we obtain 



det M(-rj) 
1 



(4(7 kf 



e 4a V _ 2e 2a V cogh 2,777 + 1 (e 2m? - l) 2 

A\\A 2 2 7, Ai 2 A 2 \ 



a 2 — a 2 



a" 



, (34) 



where the Aij denote the components of a (complex, symmetric) matrix A 
that depends on the components of the Lindblad operator as well as on the 
system parameters Ti and k. For the reader's convenience we give here the 
explicit expressions of the matrix elements, 



^22 

A 12 



= k 2 ^X>- 
= k 2 (\'- 



k\J 
U-a 

k\jt! 
H + a 



+ k 2 X" 



+ k 2 X" - 



kfj," 
H-a 

kfi" 
H + a 



k 2 (A' 2 + A" 2 + /i' 2 + // 2 ) - 2kH (AV + AV) = A21 ■ 



While 



AioA 



12^21 



[k 2 X >2 - 2kHX'fi' + k 2 fi' 2 

+ k 2 X" 2 -2kHX"fi" + k 2 fi" 2 ]- 



(35) 



is real and positive for a general non-vanishing Lindblad operator, A11A22 = 
A 12 A 2 \ + (2aka) 2 is real, but one cannot determine its sign easily for modes 
inside the horizon, where a 2 < 0. 



3.3 Modes outside the Hubble radius 

Let us first treat a mode k of the cosmological perturbations outside the 
Hubble radius (a 2 > 0), that is, before the second Hubble-radius crossing. 
Since the coupling to the environment can only lead to entanglement and not 
to disturbance in that situation, we take the dissipation coefficient a > to 
be very small, that is, a <ti \a\. We note that for such modes, 

Since we have for inflation, 
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we get arj P3 2Ti,i] e — 1 — > 1 at the end of inflation. Then, the condition 
gives, using a <C a, 

1 ~ 3^ ( cosh2 ^- 1 _ 2 
~ Pa 2 1 2a 2 ' 



yielding |?7p| ~ ^Jk/\Ai 2 \. Choosing A" = = ji" in order to implement 
a = 0, we have 

A 12 ss £; 2 (A' - /i') 2 + 2AVife(A; - W) « -2X'fi'kH , 

where in the last step we have assumed that the term with kTi dominates. 
This then gives |^ p | ~ | A'/x' | 1 and, consequently, for the corresponding pos- 
itivity time: 

t p ~ i7r 1 ln iyfl|AVI ■ (36) 



a 

At this moment, the entropy reaches the value 

iff 1 1 AV | 



S « Hit p P3 In ■ 

a 

Note that the quantity = ao/|A'^'| is invariant under arbitrary rescaling 
of spatial coordinates and has the dimension of time. 

The most striking feature in this result is its approximate independence 
of the coupling parameters (they enter t p only logarithmically). In ( fl3|) the 
influence of the environment was modelled by the parameter £, effectively 
describing localization in the field-amplitude basis as one special case of in- 
teraction. So, the result just derived supports the earlier claim in [S] that 
the details of the coupling of primordial fluctuations to the environment are 
not important. The independence of the details on the coupling is a general 
feature of systems characterized by a 'Lyapunov exponent', independent of 
whether the system is chaotic or (as is the case here) just classically unstable 
|31j . The positivity time (|36p is of the same order as the decoherence time 
discussed in [HI [TJ, [8] - Both times are generally related [17]. In this limit 
the density matrix can thus be decomposed into Gaussian pointer states 
according to ( 1271) . For strong coupling to the environment which leads to 
decoherence, these pointer states are narrow Gaussians with respect to the 
field amplitude. We note that concrete models in which a specific interac- 
tion in the action is chosen lead to results in accordance with our general 
expressions for decoherence time and entropy [H2, [50] . 

We also want to comment on the situation far outside the Hubble radius, 
but in a radiation dominated universe. There the scale factor obeys 

<V) = -rh = «e (t) , (37) 
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and thus one has again o~rj ~ 1. Repeating the above calculations now yields 
for the decoherence time, 



Hial 




(38) 



2(Ay) 2 



2 



For Hit-L ^> 1, where here t^ — o e /|A'/i'|, this is a much longer time than the 
positivity time during inflation, which means that decoherence is here much 
less efficient than during inflation, in accordance with our earlier result p]. 

The physical explanation of these results is the following. As was em- 
phasized in [3], omission of the decaying mode is sufficient to get quantum 
decoherence. However, this mode should be really small to justify this pro- 
cedure. In the super- Hubble limit, the only process making the decaying 
mode small is the expansion of the Universe; interaction with the environ- 
ment destroys the initial correlation between decaying and growing modes 
but typically makes the amplitude of the decaying mode larger (that is re- 
flected in the growth of the entropy S). During inflation, the expansion of 
the Universe is much faster (quasi-exponential) than during the radiation 
period (13T|) . Thus, the decaying mode decreases much faster during infla- 
tion, which explains the difference between ( 1361) and (1381) . Also, due to the 
exponential decrease of the decaying mode amplitude, coupling parameters 
describing real physical decoherence enter only logarithmically in t p . Note, 
however, that the fact that Ht p in the second case is larger does not mean 
that the amount of entropy gained during decoherence is less in that case, 
too. This requires special investigation. 

It is also important that Ht p > 1 in both cases. This means, in particular, 
that after the time period At ~ t p after the first Hubble radius crossing, the 
characteristic rms amplitude of the decaying mode of y is much less than the 
rms value of y in the initial vacuum state. This holds because the ratio of 
decaying to constant ('growing') mode is exp(— 3H(t— £&)), and thus this ratio 
becomes small even for a only logarithmically large exponent. Therefore, the 
Wigner ellipse remains squeezed in one direction below its vacuum width 
even after the decoherence. This is another form of the result of Sec. 2 that 



3.4 Modes inside the Hubble radius 

Let us now examine the situation inside the Hubble radius. The great gener- 
ality of f)28p makes it hard to read off the qualitative features of the positivity 
time. (See, however, the argument in [IS] that one only has to analyze one 
special example of the elliptic system in order to obtain its general qualitative 
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behaviour.) We want to consider a representative (and realistic) environment 
that is also used in [46J: a thermal bath of photons with average occupation 
number n. This calls for two Lindblad operators defined as 



1 ~ \^/a(nTY)) ' 11 -{ ) 



(39) 

h ~ 1 four 2 ~\ o 



The non-unitary part of (1281) . which is given by the dissipation coefficient 
a > 0, and the matrix A then become a sum over the different Lindblad 
operators. We then get 

A n A 22 = {(2n + l)2akH} 2 and A 12 A 21 = [(2n + l)2ak 2 ] 2 . (40) 

The determinant (I34p can now be written as 

det M(-r)) 
[(2n + l)kf 
4|a| 2 



( e 2 ^ - l) 2 (e 4 "" - 2e 2m? cos2|a|r/ + l) 

a 2 + M k 



(41) 

with |cr| 2 = k 2 — H 2 . Again assuming the coupling to be small, a <C \a\, the 
prefactor of the second term becomes much smaller than unity. If we exclude 
the unrealistic case that e 4m? — 2e 2cw? cos 2|cr|?7 + 1 ^> (e 2ari — l) 2 , the second 
term can safely be neglected to give 



det M(-ri) 



This yields a positivity time 



(2n + l)k (c2a „ ^ 2 



(42) 



1 / cr 1 \ 1 / 1 \ 
/7 P = — In l + V « — In 1 + , 43 

where the last approximation is valid for modes long after the second Hubble 
radius crossing when \a\ ~ k. (It can easily be checked that the above line 
of thought as well as the result follow from \a\ ~ k alone.) Using Friedmann 
time, the positivity time (14"3"1) reads 

a 2 . 2 ( . . 1 



( ^i«^ ln T + ^ti)- (44) 
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The positivity time is independent of the Hubble parameter Tt 7 but strongly 
depends on the (small) coupling a, in sharp contrast to the situation outside 
the Hubble radius. This is completely natural from the physical point of view 
since, inside the Hubble radius, there is no division of modes into growing 
and decaying ones: both linear independent solutions of the wave equation 
corresponding to the Hamiltonian ([T]) oscillate with amplitudes adiabatically 
decreasing with the expansion of the Universe; no more squeezing occurs. 
Thus, in this regime, the Universe expansion may not help us in getting 
decoherence, only real dissipative processes work in this direction. 

It is interesting to note that the limit n — > does not lead to t p /t c — > oo, 
but rather to t p /t e — > (In 2 2)al/16t e a 2 . This does not seem to be an artifact 
of the approximations: it is due to the vacuum energy of the photon bath, 
and may vanish after proper renormalization. 

Other concrete models of interaction may use an interaction of the form 
V = h lk (p i i(f) j k where h lk is a metric perturbation and plays the role of the 
system, while <fi is an effective ('phonon') field describing small-scale pertur- 
bations of the thermal background and plays the role of the environment. 
Clearly, such an interaction results from the kinetic term of this field. 

4 Conclusions and discussion 

Using methods from the quantum theory of open systems, we have shown 
that primordial fluctuations decohere when their wavelength becomes much 
larger than the Hubble radius and that then the pointer basis is given by nar- 
row Gaussians which approximate the field amplitude basis. Consequently, 
the entropy per mode is smaller than half the maximal entropy. Thus, the 
density matrix remains squeezed in one direction as compared to the pure 
vacuum state. However, this does not preclude the positivity of the Wigner 
function. If a particular interaction between the modes and other fields is 
given, one can calculate the Lindblad operators and all physical quantities 
in terms of this interaction and find the positivity time t p after which the 
Wigner function becomes positive everywhere and, therefore, may be approx- 
imated by a classical distribution in the phase space. In the super-Hubble 
regime, t p explicitly depends on the Hubble time if -1 and is during inflation 
only logarithmically larger than it. On the other hand, in the sub-Hubble 
regime after the second Hubble radius crossing, t p is independent of H and is 
totally determined by dissipation processes. As for the width, this leads to an 
expression analogous to (but more complicated than) the width (h/mD) 1 ^ 
for the localization of a particle. 

This crucial difference between decoherence mechanisms in the sub- and 
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super-Hubble regimes removes, we suppose, doubts regarding the possibility 
of (partial) decoherence during inflation expressed in the recent paper [19| . 
The answer to the question in this paper why we do not see similar decoher- 
ence in the Minkowski space-time is, first, because there is no super-Hubble 
regime in this case. Second, it was assumed in [19] that the 'environmental' 
(short-wavelength) modes are in their adiabatic ground state. Fields in their 
ground state are usually not able to exert a decohering influence [5J. So, the 
analysis presented here assumes that some fields are present which are not 
in their ground state. The concrete examples of such fields are cosmological 
perturbations themselves, scalar perturbations and gravitational waves, with 
scales of the order of the Hubble radius (the case not considered in [19]). But 
even if this was not true, the analysis in [3] has shown that it is impossible 
in practice to distinguish the squeezed state of the modes from a classical 
stochastic ensemble, cf. also the gedanken experiments discussed in [7| 126]. 

At last, it should be emphasized that we do not consider decoherence of 
generic long-wavelength modes but specifically only the decoherence of very 
strongly squeezed states. Then an exceedingly small interaction is already 
sufficient for a loss of quantum coherence, which has S > 1. This fact, namely 
that, even in the Minkowski space-time, strongly squeezed states are much 
more fragile than, for example, coherent states, is well known in standard 
quantum mechanics and represents the main obstacle to generate strongly 
squeezed states in the laboratory, see for example Sect. 3.3.3.1 in [3]. 

Still it should be noted that there is an agreement between our paper and 
[T9] regarding the pointer basis in the super-Hubble regime. On the other 
hand, we would not fully support the statement of [18] that decoherence is 
"extremely effective" during inflation. What follows from our results is that 
decoherence, though quick and sufficient to reach the positivity of the Wigner 
function, is not sufficiently effective, for example, to make the Wigner ellipse 
exceeding its vacuum value in all directions (the latter would correspond to 
S > Sm&x/'Z)^ Thus, our results regarding the degree of decoherence reached 
during inflation are, in some sense, intermediate between those of [18] and 
[19] . This discussion shows how subtle is the problem of decoherence and 
quantum-to-classical transition for cosmological perturbations. 

4 Note also that the ratio of inflaton gravitational to self-interaction is overestimated 
in [18) . One should take into account that, during inflation, the gravitational potential $ 
describing scalar perturbations is not constant, but slowly growing from zero at the first 
Hubble radius crossing up to its final value after the end of quasi-exponential expansion of 
the Universe. Thus, its value during inflation is much less than that after inflation which 
can be observed now. Proper account of this fact makes the inflaton gravitational and 
self- interactions of the same order. 
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